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ABSTRACT 

An  analysis  of  the  use  of  planing  data  to  predict  impact 
loads  of  a  dihedral  hydrofoil  suitable  for  seaplane  installation 
is  presented  and  compared  with  experimental  model  tests. 

The  experimental  data  verifies  the  theory  very  well  for  the 
prediction  of  maximum  loads.     No  empirical  correction  factors  are 
used  and  the  theory  agrees  with  experiment  to  an  average  error  of 
6.2$  for  the  U7  runs  analyzed.     However,   experimental  time-to- 
peak  tends  to  lag  the  theoretical,   and  the  experimental  draft 
time-history  diverges  from  the  computed  after  the  point  of  maxi- 
mum load. 

One  initial  flight-path  angle  is  considered  for  a  comparison 
between  the  maximum  loads  predicted  by  these  tests  for  a  hydrofoil 
equipped  seaplane  and  the  maximum  load  expected  of  a  chines-dry 
landing  hull. 

One  initial  flight-path  angle  is  considered  to  compare  the 
relative  efficiency  of  the  two  models   tested  as  energy  absorbing 
devices • 

In  addition,  a  sample  digital  computer  program  useful  in  the 
application  of  the  theory  to  the  prediction  of  impact  time- 
histories  is  presented. 

This  investigation,  conducted  at  the  Experimental  Towing 
Tank,  Stevens  Institute  of  Technology,  Hoboken,  N.J.,  is  the 
thesis  submitted  by  the  authors  in  partial  fulfillment  of  the 
requirements  for  the  degree  of  Master  of  Science. 


NOTATION 

M     Mass  of  impacting  body 

m     Virtual  mass  of  displaced  fluid 
w 

z  Initial  vertical  velocity 

z  Instantaneous  vertical  velocity 

z  Vertical  acceleration 

p  Mass  density  of  water,  1,938  slugs/ft. 

c  Wetted  half -width 

£  Length  of  float 

P  Wedge,  dihedral  or  deadrise  angle 

m     Virtual  mass  of  flow  plane  at  step  per  unit  distance 
in  keel  direction 

t      Geometric  angle  of  trijn 

n     An  area  of  fluid  in  a  slice  taken  perpendicular  to 
the  keel 

y      Displacement  in  a  direction  normal  to  the  keel 

y      Velocity  normal  to  the  keel 

y     Acceleration  normal  to  the  keel 

y     Displacement  of  the  step  from  origin  in  a  direction 
perpendicular  to  the  keel 

x     Displacement  in  a  direction  parallel  to  the  keel 

x     Velocity  parallel  to  the  keel 

x     Acceleration  parallel  to  the  keel 

x     Displacement  of  the  step  from  the  origin  in  a  direction 
parallel  to  the  keel 


J     Draft  in  y  direction 

Draft  of  step  in  y  direction 
F     Elemental  force  normal  to  the  keel 

ey 

F     Force  normal  to  the  keel 

y 

M     Total  mass  affected  by  passage  of  body 

Mw  =  p  cot  t   f   nw  df 
^  o 

z  Displacement  in  a  vertical  direction 

\ 
M 


"* 


F  Force  in  vertical  direction 
z 

t         2 

A^J  ^^  COS   T 

V  Horizontal  velocity 

At  An  arbitrary  increment  of  time 

^  Initial  flight-path  angle 


dt 


Variation  of  the  velocity  potential  with  time  in  the 
general  non-steady  Bernoulli  equation. 


A  Study  of  the  Hydro dynamic  Impact  Loads  on 
a  Dihedral  Hydrofoil  Suitable  for  Seaplane 

Use 

INTRODUCTION 

The  hydrofoil  as  a  seaplane  alighting  device  has  received 
serious  consideration  as  a  possible  means  of  relieving  rough- 
water  and  high  speed  landing  operational  problems. 

Compared  to  the  conventional  hull,  the  small  area  of  the 
hydrofoil  offers  a  considerable  reduction  of  landing  impact 
forces,  while  the  fixed,  relatively  short  hydrofoil  chord  limits 
center-of-pressure  travel  with  a  consequent  reduction  of  undesir- 
able pitching  moments.  To  a  large  extent  these  particular  advant- 
ages of  the  hydrofoil  are  shared  with  the  hydroski,  but  the  hy- 
drofoil offers  the  additional  advantage  of  lift/drag  ratios 
superior  to  the  ski,  especially  at  low  speeds. 

Although  the  hydrofoil  offers  much  promise  for  seaplane  use, 
with  the  exception  of  a  single  experimental  report  of  a  classi- 
fied proprietary  nature,  no  studies  of  a  general  nature  concern- 
ing its  impacting  properties  are  known. 

The  object  of  this  work  is  to  provide  a  rational  foundation 
for  the  study  and  design  application  of  impacting  hydrofoils  for 
seaplanes.  To  accomplish  this  end,  an  examination  will  be  made 
of  the  use  of  planing  data  to  compute  impact  loads  and  the 


practicability"  of  extending  its  use  to  the  specific  case  of  the 
hydrofoil;  then  to  analyze  and  compare  such  analytical  results 
with  experimental  model  tests. 

This  investigation  has  been  conducted  by  Lieutenents  Richard 
Ac  Hoffman  and  Arthur  C.  Derrick,  U0S.  Navy,  Naval  Postgraduate 
School  students  at  the  Experimental  Towing  Tank,  Stevens  Institute 
of  Technology,  Hoboken,  New  Jersey,  Special  acknowledgement  is 
due  to  Professor  B0V0  Korvin-Kroukovsky,  Mr,  Daniel  Savitsky, 
Mr,  P.  Ward  Brown  and  Mr.  Robert  Van  Dyck  of  the  Experimental 
Towing  Tank  Staff  for  their  most  helpful  interest  and  guidance  in 
this  work.   Thanks  are  also  in  order  to  all  of  the  other  members 
of  the  staff  who  so  freely  gave  of  their  assistance,  experience 
and  advice  during  this  investigation. 
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EXAMINATION  OF  THE  USE  OF  PLANING 
DATA  TO  COMPUTE  HYDROFOIL  IMPACT  LOADS 

Historical  Background 

The  earliest  impact  theories  for  the  landing  of  seaplanes 
were  based  on  the  work  of  von  Karman  (Ref.  l)  in  which  he  con- 
sidered a  wedge  shaped  body  dropped  vertically  into  a  horizontal 
water  surface.  With  the  assumption  that  momentum  was  transferred 
to  a  hypothetical  associated  mass  of  water  attached  to  the  sea- 
plane and  that  the  total  momentum  of  the  system  remained  constant, 
von  Karman  was  able  to  write  the  momentum  equation: 

Mz  =  M  z  +  m  z  (l) 

o         w 

where   M   -  mass  of  the  impacting  body 

m   -  virtual  mass  of  displaced  fluid 

z   -  initial  vertical  velocity 

z   -  instantaneous  vertical  velocity 
With  the  further  assumption  of  two  dimensional  flow  in  planes 
normal  to  the  keel,  von  Karman  considered  at  each  instant  the 
virtual  mass  is  equal  to  the  mass  of  water  contained  in  a  semi- 
cylinder  of  length  equal  to  the  length  of  the  float  and  of 
diameter  equal  to  the  width  of  the  wedge  at  the  water  surface 
(Fig.  l).  This  expression  for  the  virtual  mass  substituted 
into  equation  (l)  led  to  the  force  equation: 


k. 


z   P<rc2lcot  p 

F  =  -2 r-  (2) 

z  2 

(1  +  P^c  I) 

where  c  -  wetted  half  width 
I  -  length  of  body 

Subsequently  Wagner  (Ref.  2)  took  into  account  the  wave 
rise  at  the  wedge  surface  and  defined  the  diameter  of  the  semi- 
cylinder  of  virtual  mass  on  the  basis  of  deadrise  and  draft 
(Fig,  la);  and  Pabst  (Ref,  3)  introduced  an  aspect  ratio  factor 
to  account  for  end  losses  in  approximating  oblique  impact  on  the 
basis  of  the  two-dimensional  calculation. 

Mayo  (Ref,  U)  in  19l;£  analyzed  the  basic  von  Karman  theory 
and  concluded  it  erred  in  assuming  the  total  momentum  of  the 
system  remained  constant  in  that  it  did  not  account  for  momentum 
transfer  onto  the  wake.  Mayo  pointed  out  that  a  more  realistic 
force  equation  would  be 

Fy  =1<"»^"S^  (3) 

where  y  -  velocity  normal  to  keel 

x  -  velocity  parallel  to  keel 

m  -  virtual  mass  of  flow  plane  at  step  per  unit  dis- 
tance in  keel  direction 

in  which  the  rate  at  which  momentum  is  imparted  to  the  downwash 

is  included. 
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Mayo,  using  equation  (3),  was  able  to  apply  the  von  Karman 
theory  to  chines  dry  oblique  impact  by  imagining  the  fluid  be- 
neath the  wedge  divided  into  slices  of  length  dx  by  planes 
drawn  normal  to  the  keel.  As  the  wedge  moves  through  these 
slices,  the  flow  in  each  slice  is  assumed  to  be  identical  with 
the  flow  in  the  case  of  the  wedge  dropped  vertically.   By  assum- 
ing the  fluid  in  each  slice  to  be  independent  of  the  motion  in 
adjacent  slices,  the  total  force  on  the  wedge  is  equal  to  the 
sum  of  the  forces  of  all  the  slices  acting  on  the  wedge  plus  the 
rate  at  which  momentum  is  imparted  to  the  downwash  by  flow  planes 
sliding  off  at  the  step. 

Mayo  pointed  out  that  in  planing  the  first  term  of  equation 
(3)  is  0  and  the  hydrodynamic  force  is  equal  to  the  last  term 
and  therefore  planing  data  might  be  used  to  evaluate  the  momentum 
of  the  flow  in  the  normal  cross  sections. 

Subsequently,  in  19U8  Steiner  (Ref.  £)  followed  Mayo's  ap- 
proach to  obtain  correlation  between  planing  data  and  one  of  the 
virtual  mass  expressions,  using  limited  data  mostly  confined  to 
the  chines  dry  case.  In  1951  Smiley  (Ref.  6)  used  planing  to 
obtain  a  correlation  between  planing  and  impact  for  a  flat  plate, 
obtaining  an  impact  pressure  coefficient  in  terms  of  an  "equiva- 
lent planing  velocity". 

With  the  advent  of  heavily  loaded  high  length-beam  ratio 
hulls  which  immersed  chines  in  landing  and  of  flat  bottom  skis, 
certain  inherent  weaknesses  of  the  classic  theory  precluded  its 
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use  for  these  modern  applications.   For  example:   for  a  flat 
bottom  the  force  equation  (2)  yields  an  infinite  force j  while 
for  the  case  of  chines  wet  impact,  the  virtual  mass  becomes 
constant  at  chine  immersion  and  does  not  increase  with  increasing 
draft. 

In  19^2  Schnitzer  (Ref.  7)  treated  the  case  of  chines  wet 
impact  by  constructing  a  composite  deflected  mass  from  the  de- 
flected (virtual)  mass  apparent  at  chine  immersion  plus  a 
deflected  mass  derived  from  Bobyleff's  theory  of  the  flow  about 
a  two-dimensional  V-shape  of  finite  width.  However,  Schnitzer 's 
theory  was  lacking  somewhat  in  generality  in  that  it  attempted  to 
absolutely  define  the  geometry  of  the  virtual  mass. 

Also  in  19^2  Smiley  (Ref.  8)  presented  an  analysis  for  the 
landing  impact  in  which  the  landing  loads  and  motions  are  ex- 
pressed in  terms  of  the  planing  properties  of  the  body.  However, 
Smiley  reverted  to  an  empirico-theoretical  virtual  mass  expression 
at  chine  immersion. 

In  19!?U,  Brown  (Ref.  9)  exploited  to  the  full  the  concept 
that  impact  characteristics  may  be  predicted  with  the  use  of 
planing  data,  avoiding  the  unverifiable  definition  of  the  geome- 
try of  the  virtual  mass. 

Brown's  approach  is  especially  attractive  for  the  computa- 
tion of  impact  loads  on  such  complicated  shapes  as  hydrofoils  and 
skis,  since  it  is  unhampered  by  empirico-theoretical  expressions 
for  virtual  mass.   In  its  published  form  however,  the  derivation 
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of  Brown's  theory  breaks  down  when  the  leading  edge  of  the  body- 
becomes  wetted  as  in  an  impacting  hydrofoil.  During  the  course 
of  this  work  Mr.  Brown  became  associated  with  the  Experimental 
Towing  Tank  and  through  his  cooperation  the  derivation  was  modi- 
fied so  that  the  theory  is  applicable  to  any  impacting  body  for 
which  planing  data  is  known. 

Because  of  its  logical  simplicity  and  its  applicability  to 
any  planing  body,  Brown's  approach  is  the  one  followed  in  this 
analysis  and  the  revised  derivation  is  discussed  below. 

Seaplane  Impact  Theory  as  Applied  to  the  Hydrofoil 

The  following  analysis  considers  the  impact  of  a  straight- 
keeled  planing  body  of  arbitrary  cross  section. 
The  following  assumptions  are  made: 

1.  Trim,  %  ,  remains  constant  throughout  the  impact. 

2.  Inertia  forces  predominate. 

3.  The  body  is  a  weightless  mass,  weight  supported  by  wing 
lift. 

U.  The  virtual  mass  M  ,  is  some  function  of  the  draft. 

Consider  the  impacting  body  shown  in  Fig.  2.  In  the  fixed 

plane  AA'  a  cross  sectional  area  of  water  n   is  influenced  by 

w 

the  passage  of  the  body. 

The  momentum  of  the  fluid  in  the  ^ 

I    =  p  nw  ^s^X      ^ 
element  n   normal  to  the  keel 

w  J 
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The  area  n   is  considered  to  be  a  function  of  the  draft  f  . 
since  in  a  body  of  known  geometry  at  fixed  trim,  any  cross- 
sectional  area  may  be  determined  if  the  draft  is  known.  When  the 
body  has  passed  through  AA'  ,  the  momentum  in  the  element  is  con- 
sidered as  shed  to  the  wake.  This  concept  is  discussed  more  fully 
in  Monaghan  (Ref.  10). 

The  rate  of  change  of  momentum  ] 

\     m  Ifc  (p  n    7    Jx) 

in  an  element  fixed  in  space     J 

dn 

dn         dn  ,  ^ 

.     .  w  w       d  J 

but   St'Tf^x 

and  from  Fig.    2  J  =  y     -  x  tan  % 

5 

and  since  the  plane   is   fixed     x     is  constant 

2     ^ 
and  (£)  becomes       =  p  £x  (yg  n^  +  yg       ^-3-)  (6) 

Since  Rate  of  change  of  ^  [    Force  exerted  on 

momentum  in  an  >     =      <    body  by  element 

element  J  [ 

o  dn 

V^'^.-v*'-  af  (7) 

This  determines   the  elemental  force  acting  on  the  body  in  a 
fixed  plane.     To  determine  the  total  force   acting  on  the  body  at 
any  one   instant  of  time,   consider  the  body's  position,   velocity 
and  acceleration  to  be  instantaneously  fixed.      Now  let  the 
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previously  fixed  plane  move  along  the  body,  recalling  as  x 
varies,  J   varies  and  with  J    ,  n  . 


w 
Since  JP    =  y  -  x  tan  i  ,  under  the  above  conditions 

5 

d  J°      =  -  tan  t  dx     and 

f     S  •»  •      P     ^Tv 


X 

s 


-fs  ^  /s  dn 


=  p  ys   cot  T    \      nw  dj      +  p  ys     cot  t    ^      jf  df        M 

o 
Now  define     M    ,    the  total  mass  effected  by  the  passage  of 


the  body  as : 


M     =  p  cot  t    \       n  df 

w       r  \         w  J 


dM  U  dn 


.".  jj?  =  p  cot  t      ^  djr  (9) 

Jo 

Equation  (9)  may  not  be  readily  apparent,   but  it  may  be  proved 

by  a  series  expansion  of     n       as  a  function  of  f    • 

w  -' 

Therefore,  substituting  (9)  into  (8) 

Fy  "  ^s  Mw  +  y5     dfs  (10) 

From  Fig*    2,      J     =  z  sec  t 

d  f       =  sec  t  dz 
J    s 
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and  equation  (10)  becomes 


o  dM 

F  =  y  M  +  y   -*-=■  cos  t  (ll) 

y   •'s  w   ^s   dz 


From  Newton's  Second  Law 


F     =  -  M  y 
y  °s 


,  dM 
2       w 


.-My     ■  y     M     +y       -3 —  cos  1 
Js       ^sw      Js      dz 


o  dM  , 

••  •     y        w  /       \ 

^s  "  ^s     1?  C°S  T M"  (12) 

M(l  *  f ) 

M  dM 


M  dM 

Define     -rr  =  ^y    ,    then  -r—  =  M  ^- 

M  c    '  dz  dz 


Substituting  in  (12) 

-  7s  =  ys  T?  cos  T  xr^T  (13) 

d^t 
Determination  of   dz  and  ^ 

Consider  planing  as  a  special  case  of  impact,  Fig.  3.  In 

this  case 

y  -  0         F  =  F  cos  t 
s  z    y 

« 

y  =  V„  sin  t   where  V   is  the  horizontal  planing 


velocity 
and  substituting  in  (ll)  yields 


t-i         tt  2         2  2       w  d-H. 

F     =  VTT     sin     t  cos      t  M  -t-^- 
z         H  dz 


11. 


&2L  = 


_z_g_ 


dz        „  2      .    2  2 

VTT     sin     t  cos     t  A 
H  o 


(111) 


*.   = 


II. 


2         2  2 

Vu     sin     t  cos     t  A      J 
rl  o       o 


F     dz 
z 


(15) 


The  constant  of  integration  is  zero  since  a^    ■  0  at  z  =  0 

Thus  the  virtual  mass  ratio  --^   is  defined  in  terms  of  the 
integrated  lift/draft  curve  and  constants. 

In  order  to  change  the  equation  of  motion  (13)  which  is  in 
terms  of  velocities  and  accelerations  normal  to  the  keel  to  the 
more  useful  terms  of  vertical  velocities  and  accelerations  (as 
are  measured  in  the  towing  tank),  substitute  in  (ll) 


F  cos  t  =  F 

y      z 

y  =  V„  sin  t  +  z  cos  t 

S  n 


y  =  z  cos  t 


dM 


"2       /  •      \2   w    2 

F  =  z  cos  t  M  +  (Vu  sin  i  +   z  cos  t)  -5—  cos  t 
•  •   z  w     H  dz 

Since  F  =  -  M~    combining  with  (16)  and  (lit)  and 
z      z 

rearranging  yields 


(16) 


-  z  =  F. 


V  sin  t  +  z  cos  t 


V„  sin  t 


M(l  +  ^COS^  1) 


(17) 


12. 


2 
Letting  ^  '  =  ^  cos  t  yields 


z  =  F 
z 


V„  sin  t  +  z  cos  t 

V„  sin  t 


M(l  +-^D 


(18) 


The  resulting  non-linear  equation  (18)   does  not  lend  itself 
to  a  closed  solution,    therefore  a  step-by-step  integration  is 
necessary  to  compute  the   time  history  of  impact, 

Step-by-Step  Integration  of  Equation  of  Motion  (18) 

The  purpose  of  this  integration  is  to  determine  the  time- 
history  of  hydrofoil  impact  after  striking  the  water  under  known 
initial  conditions.     During  this  impact  the  trim  is  assumed  to 
remain  constant,   as  is  the  horizontal  velocity. 

The  horizontal  velocity,  trim  angle,   initial  vertical  veloc- 
ity, weight  of  the  body  and  a  curve  of  steady  state  planing  lift 
versus  draft  for  the  horizontal  velocity  of  impact  are  known. 
Then  to  perform  the  integration  the  only  quantity  needed  is  a 
time  interval     (At)     small  in  respect  to  the  time  from  contact  to 
maximum  load. 

If  zero  time  is  the  instant  of  contact,   the  average  draft 
during  the  first  interval  is 

e 

Z      At 
daVl  ~  — 2~ 

Then  from  a  curve  of  planing  lift  versus  draft  for  the  foil, 

the  lift  (F  in  equation  18)  for  the  average  draft  is  determined 
z 


and  from  an  integrated  plot  of  F   versus  draft,  the  quantity 

z 
-z 


£ 


F  dz  used  in  the  determination  of  u1  is  found, 
z 


To  compute  the  acceleration,  the  planing  lift  is  multiplied 
by  the  quantity 


A. 


VH 

sin 

%  + 

• 

z 

0 

COS 

T 

• 

X 

sin 

T 

1  + 


g 


2    2 

A  V„  sin  t 
o  H 


i 


F  dz 

z 


all  quantities  of  which  are  now  known. 

The  vertical  velocity  at  the  end  of  the  interval  is 


z_  =  z  -  z  At 
1    o 


and  the  draft  is 


z  +  z, 

di =  — r~ At 


Using  z..  and  d-   as  initial  conditions,  the  average  draft 
for  the  next  interval  would  be 


z1   At 
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av. 


and  the  procedure  is  repeated. 

The  most  efficient  manner  to  accomplish  the  integration  is 
in  tabular  form,   a  sample  of  which  is   given  as   Table  I.     However, 
even  at  best  the  step-by-step  integration  is  a  tedious  and  time- 
consuming   task,    and  since  each  result  is  dependent  upon  the 
previous  one,    it  is  a  process  in  which  errors   accumulate. 


m. 


An  examination  of  Table  I  reveals  that  the  step-by-step 
integration  is  in  a  form  to  be  easily  programmed  for  a  digital 
computer.  This  was  done  for  the  ELECOM  100,  and  a  sample  flow 
chart  easily  adaptable  to  any  digital  machine  is  included  as 
Table  II. 

Computation  of  Impact  by  Digital  Computer 

The  sample  flow  chart  is  so  constructed  as  to  be  useful  for 
the  determination  of  impact  time  histories  of  any  arbitrary  bodies 
for  which  planing  data  is  known. 

The  main  program  is  standard  and  need  not  be  altered  if  the 

lifting  characteristics  of  the  body  are  changed.  The  F  vs.  z 

z 

rz 

and   \   F  dz  vs.  z  functions  are  in  subroutines  which  may  be 

~>o  Z 
changed  with  comparative  ease. 

With  this  sample  program  the  only  inputs  needed  in  addition 

to  the  above  subroutines  are  V  ,  z  9    sin  t,  cos  t  ,  A  ,  g  , 

no  O 

At  ,  1  and  2.   Therefore  numerous  runs  may  be  computed  in  a  very 
short  time. 

All  theoretical  curves  except  for  the  run  of  Table  I  were 
computed  using  this  program. 

EXPERIMENTAL  INVESTIGATION 
Purpose 

While  the  planing  analogy  as  a  means  of  computing  impact 
loads  has  become  well  established  in  recent  years,  Refs.  $   through 
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9  and  11,,  all  known  experimental  data  analyzed  by  this  technique 
has  concerned  hulls  or  skis  of  low  aspect  ratio  and  a  dry  leading 
edge.  Although  intuitively  it  might  be  expected  that  the  analogy 
would  hold  for  any  planing  body  for  which  the  lift  coefficients 
are  known,  the  large  aspect  ratio  and  submerged  condition  of  the 
hydrofoil  make  it  so  different  from  the  bodies  previously  studied 
as  to  indicate  a  definite  need  for  experimental  data  to  verify 
the  applicability  of  the  use  of  planing  data  to  compute  impact 
loads.  To  this  end  the  facilities  of  the  seaplane  tank,  Experi- 
mental Towing  Tank,  Stevens  Institute  of  Technology  were  employed. 

Equipment  and  Procedure 

Two  hydrofoil  models  made  from  half-hard  brass  bar  stock 
were  used  for  the  tests.     The  first  model  consists  of  a  center- 
supported  V-shaped  hydrofoil  with  a  constant  chord  of  one  inch 
and  a  length  of  each  arm  of  five  inches.      The  dihedral  of  the 
arms   is  30  degrees  from  the  horizontal.     The  bottom  of  the   foil 
is  flat  with  sharp  leading  and  trailing  edges.     The  center  sup- 
port is  1$  inches  long  and  sharp-edged  to  prevent  excessive  drag 
from  developing   as  a  result  of  spray.     Bracing   from  the  foil  tips 
to  the  center  strut  were  added  to  insure  rigidity.     See  Fig.  U. 

The  second  hydrofoil  is  the  same  except  that  the  chord,    in- 
stead of  being  constant,   has  a  taper  ratio  of  1:2  with  a  one  inch 
chord  at  the  center  and  a  two  inch  chord  at  the  ends.     See  Fig.  U. 

The  tests  utilized  the   Tank  3    "Nutcracker"  apparatus  with  a 
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Schaevitz  type  AB-2  accelerometer  with  "damped  natural  frequency" 
on  the  order  of  90  cycles  per  second  and  approximately  0.7  criti- 
cal damping.  The  standard  Tank  3  speed  and  recording  equipment 
was  used  in  connection  with  two  Hathaway  oscillograph  recorders 
for  time-variable  measurements. 

The  readings  of  heave  and  acceleration  are  picked  up  by  and 
transmitted  by  means  of  Schaevitz  electro-mechanical  differential 
induction  transformer  units  which  give  approximately  linear  re- 
sponse of  output  voltages  to  input  mechanical  deflections. 

The  Tank  3  "Nutcracker"  apparatus  (Figs.  £,  6)  consists  of  a 
pole  free  to  move  vertically  through  two  collars  housing  low- 
friction  ball  bearings  but  limited  from  horizontal  rotation  by 
means  of  a  pivoted  arm  and  elbow  linkage.  This  apparatus  has  a 
built-in  Schaevitz  unit  to  measure  heave  of  the  pole  by  means  of 
a  string  wound  around  a  pulley  on  a  worm  gear.  A  solenoid- 
controlled  heave  lock  device  is  built  onto  the  "Nutcracker".  In 
the  tests  this  device  was  tripped  electrically  as  soon  as  the 
carriage  reached  a  predetermined  steady  speed  energizing  the 
solenoid  which  pulled  a  pin  out  of  the  pole  and  allowed  the  model 
to  drop  into  the  water.  The  "Nutcracker"  was  originally  fitted 
with  a  spring  and  cam  unloader  which  provided  a  constant  lift 
force  independent  of  the  model's  vertical  position  without  add- 
ing any  appreciable  mass  to  the  model. 

However,  after  a  number  of  trial  runs  using  the  spring  and 
can  unloader,  examination  of  the  oscillograph  tape  revealed  that 
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the  unloader  induced  such  an  extraneous  vibration  into  the  system 
as  to  make  analysis  of  the  acceleration  time-history  impossible. 
Therefore,  the  weight  of  the  model,  pole  and  appendages  was  un- 
loaded by  "negator"  type  sheave  and  line  constant  force  springs 
during  the  entire  range  of  heave,  and  downward  velocity  was 
imparted  by  using  a  ten  pound  weight  on  a  pan  attached  to  the 
pole  by  a  "C"  clamp.  The  weight  was  attached  by  nylon  cord  to 
the  fixed  (in  heave)  carriage,  Fig.  7 3   so  arranged  as  to  be 
lifted  from  the  pan  at  a  fixed  point  before  the  model  hit  the 
water.  See  Fig.  8.  Thusly  the  model  was  given  downward  momentum 
by  the  weight  and  when  the  weight  was  lifted  off  the  pan  by  the 
cording,  the  model,  being  unloaded  to  a  zero  g  condition  by  the 
"negator"  springs,  continued  downward  at  a  constant  velocity 
until  it  struck  the  water.  Examination  of  subsequent  oscillo- 
graph tapes  revealed  that  this  system,  although  hurridly  impro- 
vised due  to  the  limited  tank  time  available,  did  indeed  impart 
a  constant  vertical  velocity  after  the  weight  was  lifted  off  and 
no  discernible  vibration  to  the  accelerometer. 

The  seaplane  tank,  Tank  39   where  the  tests  were  conducted  is 
313  feet  long,  12  feet  wide  and  6  feet  deep  when  the  water  level 
is  just  even  with  the  opening  in  the  standpipe  inside  the  tank. 
A  maximum  carriage  speed  of  %0   feet  per  second  can  be  attained 
with  the  selected  speed,  amplidyne  controlled,  carriage  driving 
motor.  A  tachometer  generator,  the  armature  of  which  is  driven 
by  the  carriage  driving  motor  shaft,  supplies  electrical  information 
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to  the  field  of  the  amplidyne  system  to  provide  speed  control. 
With  a  selected  voltage,  representing  a  certain  speed,  the 
carriage  speed  can  be  maintained  to  within  0.06  feet  per  second. 

The  speed  of  the  towing  carriage  is  measured  and  recorded  by 
an  electronic  tachometer  specially  designed  and  built  by  the 
General  Electric  Company.  The  accuracy  of  this  instrument  is  in 
the  order  of  a  quarter  of  one  percent  of  the  speed  being  measured, 
in  line  with  the  degree  of  accuracy  required  in  the  drive  system. 

The  instrument  receives  its  impulses  for  speed  measurement 
from  photoelectric  scanning  discs  mounted  on  the  cable-drum  shaft. 
By  suitable  selection  of  the  number  of  scanning  holes  and  use  of 
band-spread  switches,  a  small  increment  of  speed  can  be  selected 
and  expanded  to  cover  the  full  scale  of  the  instrument  so  as  to 
insure  a  high  degree  of  accuracy  of  speed  measurement. 

As  an  auxiliary  means  of  determining  carriage  speed,  a  mech- 
anical "chronograph"  device  is  used  in  which  the  distance  travelled 
by  the  carriage  during  measured  time  intervals  is  plotted  on  mov- 
ing paper  tape. 

The  models  were  installed  in  the  overhead  towing  carriage  as 
shown  in  Figs.  5,  6,  7  and  8  with  provisions  made  to  set  the  angle 
of  trim  by  a  spirit  level  protractor.  With  the  trim  angle  set  at 
zero  degrees,  adjustments  were  made  to  insure  that  there  was  no 
angle  of  roll  or  yaw.  By  lowering  the  model  into  the  water  and 
marking  the  free  surface  intersection  on  either  side,  roll  was 
considered  zero  when  the  wetted  length  of  foil  span  on  either 
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side  of  the  centerline  was  equal  to  the  other.  Yaw  was  checked 
by  attaching  a  straight  machined  bar  with  its  lengthwise  axis  in 
the  direction  of  motion  to  the  model  strut.  Then,  if  by  moving 
the  carriage  forward  the  edge  of  the  bar  remained  at  a  constant 
distance  from  the  tank  wall,  yaw  was  considered  zero. 

Calibration  in  heave  was  made  by  moving  the  pole  up  and  down 
in  one  inch  increments.  Heave  tape  readings  were  made  at  each 
position.   Therefore  inches  of  tape  per  inch  of  actual  heave 
could  be  determined. 

Calibration  in  acceleration  was  made  by  detaching  the  accel- 
erometer  from  the  pole  and  standing  it  upright  (zero  "g"),  on  its 
side  (minus  one  "g")  and  upside-down  (minus  two  "g"),  and  taking 
tape  readings  at  each  position.  Thus  inches  tape  per  "g"  could 
be  determined. 

Calibration  in  heave  and  acceleration  was  made  daily  before 
starting  runs  and  whenever  a  change  in  set-up  was  made. 

In  addition,  before  each  run  heave  was  marked  on  each  tape 
at  the  bottom  stop;  where  the  trailing  edge  tips  of  the  model 
just  touched  the  water;  and  where  the  trailing  edge  apex  of  the 
model  just  touched  the  water.   From  the  intersection  of  the  apex 
touching  trace  with  the  running  heave  trace,  the  beginning  and 
end  of  impact  was  determined. 

Runs  were  made  at  a  constant  forward  speed  of  about 
30.7  ft. /sec.  and  for  each  trim  angle,  various  vertical  veloci- 
ties, see  Table  III.   The  vertical  velocities  in  each  case  could 
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not  be  precomputed  due  to  friction  in  the  pole  and  the  stretching 
of  the  nylon  cord  which  varied  the  length  of  drop  over  which  the 
weight  acted.  All  vertical  velocities  were  computed  from  the 
heave  vs.  time  oscillograph  tape. 

During  the  planning  phase  of  these  tests  it  was  concluded 
that  in  full  scale  landings  full  ventilation  of  the  upper  foil 
surface  would  undoubtedly  occur.   In  order  to  simulate  actual  ex- 
pected conditions,  the  experiments  were  planned  to  achieve  separa- 
tion at  the  leading  edge  and  consequent  ventilation  by  a  combina- 
tion of  high  forward  speed  and  a  sharp  leading  edge. 

As  demonstrated  by  Wu  (Ref.  12)  with  full  ventilation  the 

lift  coefficient  of  a  submerged  planing  flat  plate  may  be  pre- 

2  Tf   sin  t 

dieted  by  the  Rayleigh  formula  CT  *  i  ''   i  . ,  which  can 

be  modified  for  the  effect  of  dihedral  on  angle  of  attack  to 

yield  CT  =  ■, — y^s  .n cos  K .  After  a  few  trial  impact  runs  it 

L  k   +7ysin  t  cos  p  ^ 

became  evident  that  the  modified  Rayleigh  lift  coefficient  was 
not  being  realized  at  low  trim  angles.   It  was  concluded  that 
this  was  due  to  the  fact  that  the  models,  being  designed  for 
rigidity  and  economy  of  fabrication,  had  too  large  a  leading  edge 
angle  to  achieve  separation  at  low  trims,  and  that  at  trims  below 
the  critical  separation  trim  the  upper  surface  was  actually  gen- 
erating negative  lift.  Since  funds  were  not  available  to  recon- 
struct the  models,  the  trim  angle  was  increased  to  a  minimum  of 
12°  for  the  tapered-chord  foil  (Foil  #l)  and  10°  for  the  constant- 
chord  foil  (Foil  #2),  where  positive  lift  was  attained.  It  was 
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further  concluded  separation  and  subsequent  ventilation  occurred 
at  and  above  these  trims  and  they  were  the  minimum  examined  in 
the  impact  condition. 

In  order  to  obtain  accurate  lift  data  for  the  impact  condi- 
tions, planing  data  was  obtained  experimentally  for  both  foils 
at  the  trim  angles  and  horizontal  speed  considered  in  impact. 
This  experimental  lift  vs.  draft  data  is  plotted  in  Figs.  9  and 
10  for  Foils  1  and  2,  respectively.  This  planing  data  was  ob- 
tained by  the  same  apparatus.  After  unloading  the  model  and  pole 
with  "negator"  springs  to  a  zero  "g"  condition,  known  weights 
were  attached  to  the  model  and  it  was  run  down  the  tank  at  the 
forward  speed  to  be  examined  in  impact  in  a  f ree-to-heave  condi- 
tion. When  the  predetermined  steady  speed  was  reached,  the  heave 
(model  draft)  was  read  from  the  visual  scale  of  the  oscillograph. 
The  plots  of  lift  vs.  draft  were  obtained  by  varying  the  load. 

COMPARISON  OF  THEORETICAL  RESULTS 
WITH  EXPERIMENTAL  DATA 

General  Discussion   In  order  to  determine  the  applicability  and 
degree  of  accuracy  of  the  developed  impact  theory  to  the  specific 
case  of  the  dihedral  hydrofoil,  a  comparison  was  made  with  ex- 
perimental model  test  data.  The  results  of  the  tests  are  pre- 
sented as:  (l)  experimental  curves  of  vertical  steady-state 
planing  force  vs.  draft,  Figs.  9  and  10  for  foils  1  and  2  respec- 
tively, made  at  the  single  horizontal  velocity  and  at  the  various 
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trim  angles  examined  in  impact.  The  trim  angles  used  were  high 

enough  to  insure  ventilation  in  the  planing  condition;  (2)  curves 

of  the  function   IF  dz  vs.  draft  obtained  from  graphical  inte- 

J  o 
gration  of  the  lift/draft  curves.  These  curves  are  presented  as 

Figs.  11  and  12  for  foils  1  and  2  respectively;  (3)  plots  of  maxi- 
mum loading  per  run  vs.  initial  flight-path  angle,  Figs.  13  and 
III.  for  foils  1  and  2  respectively;  (U)  time-history  comparisons 
with  results  obtained  from  the  stepwise  integration  of  the  equa- 
tion of  motion,  equation  (18).  See  Figs.  15>  through  32;  and  (£) 
loop  curves  to  compare  energy  absorbing  efficiency  of  the  two 
models  at  one  initial  condition,  see  Fig.  33« 

Experimental  Data   The  experimental  data  obtained  in  the  tests 
were:  (l)  curves  of  steady-state  planing  force  vs.  draft  necess- 
ary in  the  computation  of  the  theoretical  time-histories  and 
(2)  experimental  time -histories  of  draft  and  vertical  accelera- 
tion. Experimental  draft  and  acceleration  were  measured  directly 
from  the  oscillograph  tape.  The  initial  vertical  velocity,  meas- 
urement of  which  is  extremely  important  for  an  accurate  comparison 
between  theory  and  experiment,  was  computed  from  the  slope  of  the 
draft  trace  just  before  water  contact.  Correct  interpretation  of 
the  slope  of  the  draft  trace  is  quite  dependent  upon  the  judgment 
of  the  observer.   For  this  reason,  the  initial  vertical  velocities 
were  computed  independently  by  both  investigators  and  the  recon- 
ciled values  were  used  in  this  report. 
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Comparison  with  Theory       Time-histories  of  experimental  draft  and 
acceleration  for  18  of  the  runs  listed  in  Table  III  are  compared 
in  Figs.   15>  through  32  with  those  computed  from  equation  (l8)  for 
the  same  initial  conditions •     The  runs  compared  were  chosen  at 
random  and  are  considered  representative  of  the  three  trims  used 
and  approximately  cover  the  range  of  initial  flight-path  angle, 

v    ,    used  in  the  tests .     Again,   it  is  emphasized  that  full 
ventilation  is  expected  in  full  scale  landings  and  that  the  test 
impact  data  is  for  the  case  of  upper  surface  ventilation. 

Results  and  Discussion 

Acceleration       Examination  of  Figs.  V~>  through  32  reveals  that 
the  presented  theory  is  in  excellent  agreement  with  the  experi- 
mental values  of  maximum  acceleration.     However,   a  general  tend- 
ency of  the  experimental  time-to-peak  to  lag  the  computed  is 
apparent.     Although  this  lag  could  be  blamed  on  the  difficulty 
in  establishing  an  accurate  reference  line  for  acceleration 
measurement,   it  is  felt  that  this   time  lag  is  a  manifestitation 
of  the  inherent  weakness  of  the  planing  analogy,   i.e.,   the  treat- 
ment of  impact  as  a  quasi=static  phenomenon.      It  is  well  known 
that  pressures  on  submerged  objects   as  determined  by  the  general- 
ized Bernoulli  equation; 

contain  a     — j£     term  which  increases  in  importance  with  increasing 
o  t 
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rapidity  of  time  change.     For  example,    in  the   two-dimensional 

wedge  impact  treated  by  "expanding  plate  theory"  of  von  Karman 

and  Wagner,    the  major  part  of  the  force  results  from  the     £*■ 

d  t 

term.  However,  this  theory,  as  most  potential  impact  theories, 

treats  the  particular  case  of  uniform  penetration  velocity  in 

which  case  the  ^  term  furnishes  the  major  resisting  force. 
a  t 

In  decelerated  penetration,  as  in  a  seaplane  landing,  the  -rr 

o  t 

term  would  be  reduced  by  the  inertia  of  the  decelerating  virtual 
water  mass.  In  the  case  considered  in  the  present  report,  the 
consideration  of  impact  as  quasi-static  can  be  considered  as  hav- 
ing provided  a  satisfactory  first  order  answer.  With  increasing 
impact  velocity  however,  the  effect  of  the  -r£  term  as  affected 
by  deceleration  can  be  expected  to  become  more  and  more  important, 
manifesting  itself  in  a  lag  of  the  resisting  force.  As  penetra- 
tion velocity  is  reduced  and  approaches  zero  the  conditions 
approach  static  planing  and  satisfactory  agreement  of  maximum 
acceleration  is  expected  and  in  fact  obtained. 

It  is  felt  a  more  precise  and  rigorous  impact  theory  for 
ventilated  hydrofoils  can  be  obtained  from  potential  flow  theory. 
Green's  paper,  Ref.  13,  provides  potential  equations  for  the 
discontinuous  flow  about  a  flat  plate  as  does  the  work  of 
Rayleigh,  Ref.  lU,  and  Milne-Thomson,  Ref.  V~>.     However,  all  of 
these  references  express  the  potential  function  in  terms  of  an 
intermediate  mapping  plane  and  the  mathematical  difficulties  in 
expressing  the  potential  in  terms  of  the  physical  plane  precluded 
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further  examination  of  this  approach  in  the  time  available. 
It  is  felt  that  further  examination  of  the  application  of  the 
time  dependent  potential  flow  theory  to  hydrofoil  impact  will  be 
most  fruitful o 

Draft  Comparison  of  the  computed  and  experimental  draft  curves 
of  FigSo  15  through  32,  show  good  agreement  during  initial  pene- 
tration. However,  all  experimental  values  of  maximum  draft  are 
larger  than  those  computed.  It  is  felt  that  part  of  the  reason 
for  this  deeper  penetration  is  due  to  the  time-lag  effects  dis- 
cussed above.  Due  to  the  time-lag,  the  energy  is  absorbed  more 
slowly  than  predicted,  resulting  in  a  deeper  penetration. 

Another  reason  for  the  variation  is  considered  to  be  due  to 
separation  and  ventilation  effects.  As  the  vertical  velocity 
becomes  small  near  maximum  draft,  the  flow  tends  to  close  over 
the  upper  surface  of  the  foil,  and  for  these  particular  foils 
this  reduces  the  lift.  With  less  lift  to  resist  penetration, 
the  draft  is  increased. 

The  ventilation  effects  are  especially  evident  at  low  trim 
angles.  For  foil  1  at  12°  trim  and  for  foil  2  at  10°  trim, 
Figs.  V?   through  17  and  Figs.  2\\   through  26  respectively  show 
experimental  draft  values  diverging  from  the  computed  after 
maximum  draft,  while  at  higher  trims  the  experimental  curves 
generally  parallel  the  computed  after  maximum  draft.  Since  loss 
of  proper  ventilation  is  more  critical  at  the  lower  trims,  the 
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divergence  of  the  experimental  draft  from  the  computed  low  trim 
indicates  definite  closing  in  of  the  flow  at  low  trim  and  probable 
closing  in  effects  at  the  higher  trims. 

It  is  felt  that  these  pronounced  ventilation  effects  are  due 
in  a  large  part  to  the  poor  hydrodynamic  design  of  the  model  foils 
tested.  The  authors  of  this  work  have  proposed  tests  of  a  dihedral 
hydrofoil  constructed  for  a  dynamic  model  of  an  actual  aircraft 
to  the  Experimental  Towing  Tank,  and  the  results  obtained  therefrom 
should  be  in  much  better  agreement  with  theory. 

Variation  of  Maximum  Loading  with  Initial  Flight-Path  Angle 
Figs.  13  and  lh  reveal  that  the  maximum  loading  of  a  dihedral 
hydrofoil  increases  linearly  with  initial  flight-path  angle,  and 
with  a  slight  increase  of  maximum  loading  with  increasing  trim. 

Effect  of  the  Virtual  Mass  Term  (  1  +^ ' )   of  Equation  (18) 
Analysis  of  ten  hand  calculated  step-by-step  integrations  of  the 
equation  of  motion,  of  which  Table  I  is  a  sample,  reveals  that 
the  =-— — r  term  of  equation  (18)  is  conservative,  with  a  mini- 
mum  observed  value  at  maximum  loading  of  .90  .   Therefore  its 
omission  in  initial  estimate  of  maximum  loading,  which  greatly 
simplifies  hand  calculation,  should  result  in  a  maximum  over- 
estimate of  maximum  loading  of  approximately  10$. 

Usefulness  of  Sample  Computer  Program  of  Table  II   Comparison  of 
results  obtained  from  the  digital  computer  program  of  Table  II 
with  hand  calculated  results  revealed  the  use  of  the  computer 
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resulted  in  a  much  greater  degree  of  accuracy  than  possible  with 
a  slide  rule  at  a  great  saving  in  time.  Calculation  of  the  com- 
plete time-history  of  impact  by  hand  took  approximately  two  hours. 
Calculation  by  the  digital  computer  took  about  3  minutes  of  com- 
puter running  time. 

Possible  Application   Although  the  hydrofoil  models  of  these 
tests  were  designed  primarily  for  verification  of  the  use  of 
planing  data  to  compute  hydrofoil  impact  loads,  the  general  con- 
figuration of  the  models  is  envisioned  as  being  suitable  for 
seaplane  installation.  To  afford  some  first-order  degree  of 
approximation  as  to  the  advantages  of  such  a  hydrofoil  system, 
consider  the  models  to  be  scaled  up  by  a  factor  of  11.9.  This 
corresponds  to  an  aircraft  having  a  landing  speed  of  63  knots 
and  a  gross  weight  of  90^0  lbs.  for  foil  1  and  8000  lbs.  for 
foil  2.  Since  these  weights  and  speed  are  comparable  to  the 
Grumman  JRF,  the  £  foot  beam  and  30°  deadrise  angle  of  the  JRF 
were  selected  as  further  comparison  parameters.  Using  Milwitsky's 
tables  for  chines-dry  impact,  Ref.  16,  and  the  10.3  flight-path 
angle  of  runs  1  and  37  for  foils  1  and  2  respectively,  maximum 
load  factors  for  the  chines-dry  hull  impact  were  computed  as 
7.3"g"  and  7»6ug"  respectively.   Comparison  of  these  predicted 
loads  with  the  hydrofoil  loads  as  listed  in  Table  III  for  the 
same  initial  conditions  show  that  foil  1  provides  a  h5%   and 
foil  2  a  $6%   reduction  in  maximum  impact  load  over  the  chines-dry 
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hull.  This  reduction  is  accomplished  by  a  full  scale  foil  of 
8.£  ft.  span  and  the  following  chord  dimensions: 

For  foil  1  -  an  apex  chord  of  11.9  inches  and  a  tip  chord  of 
23.8  inches. 

For  foil  2  -  a  constant  chord  of  11.9  inches. 

To  consider  the  foils  as  energy  absorbing  devices,  the  loop 
charts  of  Fig.  33  were  constructed.  Graphically  integrating  the 
area  under  the  curves  to  the  point  of  maximum  penetration  checked 
the  initial  kinetic  energy,  -  M  z   ,  within  $%.     Furthermore, 
comparing  the  area  under  the  curves  with  the  square  area  of  an 
ideal  shock  absorber  showed  both  foils  to  have  an  ideal  effici- 
ency of  approximately  %6%.     However,  taking  the  ratio  of  rebound 
energy  (area  under  lower  limb  of  loop  curve)  to  input  energy  (area 
under  upper  limb),  showed  that  foil  1  returned  only  17«5£of  the 
original  kinetic  energy  to  the  system  as  rebound  while  foil  2 
returned  28.1$.  Therefore,  a  tapered  foil  corresponding  to  the 
dimensions  of  foil  1  would  provide  a  more  damped  landing  device 
with  lesser  likelihood  of  the  heavy  damage  experienced  with  the 
rebound  impacts  of  the  conventional  hull. 

Conclusions   An  analysis  of  the  experimental  data  obtained  dur- 
ing hydrodynamic  impact  of  a  dihedral  hydrofoil  resulted  in  the 
following  conclusions: 

(l)  In  general,  the  use  of  the  quasi-static  planing  analogy 
to  compute  hydrofoil  impact  loads  gives  a  good  all-around  approach. 
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In  particular,   computed  maximum  loadings   are  in  very  good  agree- 
ment with  experiment,   showing  an  average  error  of  6.2$  for  U7 
runs, 

(2)  Experimental  time-to-peak  tends  to  lag  the  computed. 

(3)  Correspondence  between  computed  and  experimental  draft 
is  good  until  maximum  force,   after  which  a  divergence  is  noted 
due  to  the  time-lag  of  (2)   above  and  separation  effects. 

(U)  The  maximum  impact  loading  of  a  dihedral  hydrofoil 
increases  linearly  with  initial  flight-path  angle  and  increases 
slightly  with  increasing  trim. 

(£)  The  effect  of  the  Virtual  Mass  Term  (=-| — j-1  of  equa- 
tion (18)  is  conservative,   with  a  minimum  observed  value  at 
maximum  loading  of  ,90. 

(6)  The  sample  digital  computer  program  of  Table  II  improves 
the  accuracy  of  the  calculation  of  draft  and  acceleration  time- 
histories  and  reduces  the  time  required  for  such  calculation 
from  2  hours  to  3  minutes. 

(7)  Compared  with  chines-dry  impact  at  a  10,3°  flight-path 
angle,   the  dihedral  hydrofoil  system  offers   a  reduction  of  maxi- 
mum impact  loading  in  the  order  of  $0%, 

(8)  At  approximately  8.9     initial  flight-path  angle,   both 
foils  tested  had  an  ideal  shock  absorber  efficiency  of  56%;   how- 
ever the  larger  tapered  foil  1  returned  but  17.5$  of  the  original 
kinetic  energy  to  the  system  as  rebound  energy,  while  foil  2  re- 
turned 28 .h%,   and  therefore  a  foil  similar  to  foil  1  is   recom- 
mended for  its  damping  qualities. 
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Sample   Flow  Chart  for  Step-by-Step  Computation  of  Equation 
of  Motion  by  use  of  a  Digital  Computer 
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Cell  Fill  and  Working  Storage  for  Sample  Digital  Computer  Program 
to  Determine  Time  History  of  Impact 
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